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Abstract 

Stimulating a system with time dependent sources can enhance instabihties, thus 
increasing the critical temperature at which the system transitions to interesting low- 
temperature phases such as superconductivity or superfluidity. After reviewing this 
phenomenon in non-equilibrium BCS theory (and its marginal fermi liquid general- 
ization) we analyze the effect in holographic superconductors. We exhibit a simple 
regime in which the transition temperature increases parametrically as we increase the 
frequency of the time-dependent source. 
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1 Introduction 

Considering time-dependent couplings among quantum fields leads to a rich theory and phe- 
nomenology with a number of important applications]^ One interesting effect was explained 
early on in following the development [2] of a microscopic theory of superconductivity 
as pairing of charged fermions leading to spontaneous breakdown of electromagnetic gauge 
symmetry. The work [1] argued that time dependent high-frequency electromagnetic fields 
could modify the density of states of the fermions in such a way as to increase the transi- 
tion temperature Tc of the material (without appreciably heating it up), and the work [3] 
analyzed the net effect taking into account phonon interactions. This effect had been seen ex- 
perimentally, and has been tested more recently with positive results reported [Ij. Recently 
there appeared work [5] analyzing the closely analogous effect in superfluids built from cold 
atoms, which may have the advantage of additional experimental control. Other interesting 
low-temperature phases can also be stimulated in a similar manner [6]. In each of these 
works, the calculations were controlled by treating the externally applied time-dependent 
sources perturbatively. 

Abstractly in quantum field theory, it is easy to see how time-dependent couplings can 

^not least among them the origin of structure in the universe. 
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enhance an instability such as a superconducting transition. Consider a scalar field with a 
time-dependent coupling g{t) 



S = j d^X-^^^ I (50)^ - ^"^V^ + (^interaction^ ■ 

Canonically normalizing the scalar field by setting = 4>/g{t) yields an action 



S= d^x{ {d4>y 



2 n ( '9 
m — 2 [ - \ H — 



9/ 9 



^interaction { • 



The quantity in square brackets behaves as an effective mass squared, miff. The shift is 
caused by the time-dependent coupling. For a range of functions 9{t), the shift is negative. 
In the context of a Ginzburg-Landau model, this reduction in rrigg enhances the instability, 
raising the transition temperature below which the scalar field condenses. (Depending on 
the implementation, it may also thermally excite the system, so determining the net effect 
requires more detailed analysis.) 

In practical applications, one does not a priori control g{t) in the low-energy effective 
theory; one may either control the chemical potential as a function of time, or irradiate 
a sample electromagnetically, or both. We are therefore motivated to analyze the effects 
of such time-dependent perturbations in more general quantum field theories, in particular 
strongly coupled field theories which exhibit superconducting transitions [7J. 

In §§3 and 4 we will derive simple (but controlled) analytic estimates exhibiting our 
effect, and check them with appropriate numerical solutions. By working in a simple probe 
limit and estimating the particle production and decays, we will check that for a reasonable 
class of examples the main effect of the time dependence is to raise the effective of the 
system (as opposed to heating it up)j^ In modern pump-probe experiments, one can study 
non-equilibrium dynamics, but one typically heats up the system. It would be interesting to 
model this physics holographically as well. 

Before turning to that, however, we will briefly review the effect in BCS theory and 
comment on its generalization to strongly coupled analogues along the lines of |9] . |^ 



2 Time Dependence and Pairing 

In this section we briefly review the mechanism of pQ, originally applied to BCS pairing in 
Fermi liquid theory, and comment briefly on its generalization to the non-Fermi liquid model 

^ Other interesting works on time-dependent finite density holographic field theory have appeared recently, 
including [8] . These have typically focused on a return to equilibrium, whereas in this work we are concerned 
with the effect of nonequilibrium dynamics in enhancing phase transitions. 

■^Although no complete model of high-temperature superconductivity has been derived, it is interesting 
to consider the corresponding question in those materials. Ultimately a factor of 2 in Tc would make the 
difference between current high- Tc superconductors and room-temperature ones. 
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of [9]. In BCS superconductors [2J, one calculates the effective potential for a Cooper pair; 
the resulting equation of motion is the gap equation, essentially the equation of motion for 
the charged scalar field condensate A, which takes the form 

'^'^^°'dE^^^^^^il-2n,n). (2.1) 

Here A^(0) is the density of states at the Fermi surface {E — )■ 0), A the strength of the 
pairing interaction, and npD = l/{e^^'^ + 1) the Fermi-Dirac distribution as a function of 
temperature T. For T larger than a critical temperature Tc ~ EcutoS^^'^^'^'^^^K the only 



solution to (2.1) is A = 0. Below there is a nontrivial solution for A, leading to the Higgs 



mechanism and superconductivity (see [TU] for a pedagogical review) 



The kinematics of the Fermi surface is crucial in (2.1); the effect dies away exponentially 
as the density of states A^(0) at the Fermi surface goes to zero. The finite-temperature distri- 
bution npD smooths out the Fermi surface and eventually eliminates the superconductivity. 
Thermal occupation of quasiparticle states just above the original Fermi surface prevents 



them from contributing to (2.1). Eliashberg [T] noted that time dependent perturbations 
could counteract this effect of finite temperature, by exciting the quasiparticles just above 
the original Fermi surface to significantly higher levels. This effectively shifts the distribution 



function npr, in (2.1) to a smaller value, 

npD ^ npD + with Sn < 0. (2.2) 

With the excited quasiparticles out of the way, the system can access the states near the 
original Fermi surface, enhancing pairing and condensation. Chang and Scalapino [3] and 
others studied the further effects of phonon-quasiparticle interactions, finding that additional 
enhancements are possible. 

In the next section, we analyze the effects of time dependent perturbations on holographic 
superconductors]^ In these systems, we will find that time dependent perturbations also en- 
hance Tc, even allowing for superconductivity when the average chemical potential vanishes. 
Although it goes in the same direction, the mechanism behind this is not the same in detail 
as that in the BCS example. One question raised in the earlier works (see e.g. [5]) is how to 
extend the calculations beyond a perturbative treatment in 6n (equivalently, perturbative in 
the original time dependent source) to see how large the enhancement of could become. 
In our case, the holographic methods allow us to analyze large 6n, though in an unrealistic 
regime of large-N gauge theory. As we will see, the effect can be parametrically large in this 
regime. 

An example that may be closer to the BCS case is to consider pairing and time dependent 
perturbations in the free sector of the model [5], which produces a marginal Fermi liquid 
upon mixing with a putative locally critical strongly coupled sector of fermion operators. 
Time dependent perturbations could again be used to excite thermally distributed states 
away from the original Fermi surface, increasing Tc from its time-independent value. 



''More precisely, the gravity side is dual to a superfluid, but this distinction will not be important in our 
case. Some of the nonequilibrium effects are in fact easier to observe in superfluids. 
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3 Holographic Time-Dependent Superconductors 



In this section, we analyze the effect of a time- dependent chemical potential on in a 
holographic superconductor of the type introduced in [7j|^ We will take the temperature 
T ^ Tc, a limit which is sufficient for obtaining the critical temperature. This corresponds 
to the onset of the instability and reduces the nonlinear system to a single linear differential 
equation. In ^we will study the condensed solution that arises when T < Tc. 

The basic framework is 4D AdS gravity with a U{1) gauge field and a charged scalar: 

S = j d^x^ {^+J2- I^M-^''' - 1^*^ - - ^'l-^l') , (3.1) 

where L is of order the AdS radius. The background metric corresponds to an AdS- 
Schwarzschild black hole, 

ds' = -f{r)dt' + ^ + -y^idxl + dxl) , fir) = 77 U - ^ , (3-2) 



where the temperature is 



„ 3r I 



47rL2 



(3.3) 



The dual QFT has an operator O dual to $, and its expectation value is the order 
parameter for superconductivity. The dimension A oi O and the mass of $ are related by 

A(A - 3) = m^L^ . (3.4) 

The U{1) bulk gauge field is dual to a global symmetry current in the QFT, so strictly 



speaking (3.1) describes a superfluid. This distinction, however, will not play an important 
role here, and it is possible to make the boundary gauge field dynamical. We will comment 
more on this below. 

We will consider a time-dependent chemical potential in the dual QFT: 

At -> fi(t) , for r 00, (3.5) 

and we want to find its effects on the charged scalar field $. 

A time-independent chemical potential has been analyzed in [7] and other works. The 
coupling of the charged scalar $ to the gauge potential At leads to a negative contribution to 
the effective mass squared for $(r), causing it to condense for r < r^, where Tc is the radial 
position below which the effective mass squared is sufficiently negative. In the present work 
we wish to compute the contribution of time dependent effects to the effective mass squared, 
and determine the new Vc below which $ condenses. This provides a measure of the shift of 
Tc in the dual QFT. 

^Effects of a spatially dependent chemical potential on holographic superconductors were studied in 
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3.1 Probe Limit 

To illustrate the effect, it is sufficient to work in the probe limit where and $ do not 
backreact on the metric. This entails adding stress-energy sources in the bulk which are 
subdominant to M|,/L^, where Mp is the bulk four-dimensional Planck mass scalej^ In this 
regime, the equations of motion are 

V^F^^ = iqi^d""^ - ^d"^) + 2q'^\^\'^A'' 

D^D^"^ = 771^$ . (3.6) 

To begin, consider an ansatz where the only nonzero gauge field is At. In the static case, 
the phase of $ is required to be constant, something which is no longer true once time- 
dependence is allowed; for instance. Maxwell's equation along the r-direction now reads 

dtdrAt = iqf{r){^dr^ -^dr^) . (3.7) 

The dynamics is then easier to analyze in unitary gauge, where a gauge transformation is 
used to set the phase of $ to zero. The gauge transformation achieving this is 

K = A,- ^9,(arg($)) . (3.8) 

Henceforth we work in terms of A' and drop the primes; in this gauge, we need to allow for 
both At and Aj. components. 

In unitary gauge the equations of motion simplify further. Maxwell's equations along the 
time and radial directions now read 

\dr{r'Frt) - %^'At = (3.9) 
dtFrt = 2q^f^'Ar. (3.10) 

The second equation shows explicitly that the source for Ar is given by the time and 
radial dependence of At. Note that because of our unitary gauge choice, $ is real here. Its 
equation of motion is 

From this equation, we read off an effective mass squared 

ml, = m'-^^ + q'fAl. (3.12) 



^More generally, other scales may arise in the UV completion of the system, for which similar remarks 
hold. 
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We will be interested in determining at what scale this becomes sufficiently negative to 
produce an unstable mode. 

Finally, we need to impose boundary conditions. In the static case, the boundary condi- 
tions at infinity are 

$ = ^ + ... , A, = fi-^ + ... (3.13) 

where 

A± = ^ ± ^V9 + 4m2L2. (3.14) 

The fall-off for $ corresponds to the expectation value for the operator dual to the scalar 
field, 

{O) = (3.15) 

when there is no external source/] In order to maintain a finite current squared J^J^ at the 
horizon, one sets At = there [7]. 

To generalize this to the time-dependent case, the first step is to replace /x — )■ fi{t) and 
keep similar fall-offs at infinity. In order to determine the critical temperature, we need only 
focus on solutions at or near $ = 0; we will turn to this next. 



3.2 Estimate in a Simple Tractable Regime 



We are interested in the radial scale at which the instability of the scalar field $ kicks in; 
this will determine the critical temperature Tc for the superconducting transition in the dual 
field theory. In order to estimate this, we can simply study the unstable solutions at $ = 
and determine at what radial scale m^^ goes sufficiently negative. Setting $ = in (3.9), 
(3.10), and (3.11), the solution simplifies to 

Ftr '- 



where our choice of notation will become clear momentarily. Now let us impose that 



(3.16) 



(3.17) 



as we approach the boundary r — )■ oo. 
Plugging (3.17) into (3.16) yields 

Ftr = dtAr - dr-At 



dtAr- 



p{t) 



(3.18) 



valid near the boundary. The boundary condition for At near the horizon will determine 
the relation between the chemical potential and charge density, leading to p{t) = r^^(t). 
Anticipating this, we obtain 



(3.19) 



"^In the range 



< m < m 
consider the boundary condition $ — !■ (f)^/r'^ 



''BF 



BF 



L ^ the mode is also normaUzable and we could instead 
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where we set the integration constant to zero so that on average {Ar)t = 0, as in the static 
case. For a sinusoidally varying chemical potential 



n{t) = /io — 6fism{ut) 



(3.20) 



we have 



— — COS[Ujt) 



(3.21) 



Notice that the amplitude of A,, is controlled by 6fi/uj. 

In fact, (3.17) and (3.19) solve the equations of motion everywhere in the bulk, not 



just near the boundary. The boundary conditions near the horizon determine the relation 
between /i(t) and p{t). Let us impose that — > at the horizon (so that J^J'^ is finite 
there, where is the bulk gauge current |7j). This is satisfied if we set fi(t) = p(t)/r^, i.e. 



At = pit) (^1 



(3.22) 



Now in the near-boundary regime r ^ L, the emblackening factor /(r) — )■ r'^/L'^. From 



the expression (3.12) for mlg, the charge density will enhance the instability of $ provided 



that 



^4 

Aj->Al 



(3.23) 



Because of the factor of w in (3.21), for sufficiently large u this will be satisfied, and the 



effect is enhanced by increasing the time dependence (increasing oj). 

Note that m^^ — )■ both near the boundary and at the horizon. An instability will ap- 
pear if At can make the effective mass squared sufficiently negative over a large enough range 
of radial scales. In particular, if we begin with a bare m? which is right at the Breitenlohner- 
Freedman bound, the net negative contributions to m^g of of the time dependence must 
induce a condensate. 

In §3.3| and ^ we will further analyze this, and solve the equation of motion for $ 



numerically near Tc. Its basic features can be understood qualitatively. The instability will 
set in at some which is greater than its value in the absence of the oscillating charge 
density. This increases Tc since it increases the radial position that needs to be cloaked by 
the horizon if we wish to heat the system up above the superconducting transition. These 
effects are illustrated in Figure [l| where mgfj(r) is plotted for different values of the frequency. 
Moreover, in the case /xq = we see that the time-averaged chemical potential is zero, and 
the superconductivity follows entirely from the oscillating contribution to the time average 
{9''Al)t. 



3.3 Determination of Tc 

For the purpose of finding the dependence of Tc on 5p and u it is sufficient to restrict to 
T ^ Tc. In this approximation, $ should be everywhere small and its backreaction on the 
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Figure 1: Effective mass mgg(r) for different values of the frequency w. In the plot, m'^L'^ = —2, 
r^/L = 0.5, and for reference the BF bound m?^pLp' = —9/4 is shown. Increasing the frequency 
makes m^g more negative over a larger range of radial coordinate. 



gauge fields (which is quadratic in $) may be ignored. and Ai are then given by (3.19) 



and (3.22). 



Furthermore, we focus on the large frequency regime, and approximate the time-dependent 
equation of motion for $ by its time average. This reduces to the ordinary linear differential 
equation 



where the effective mass (3.12) becomes 

,2 



(r_ 

f 



[^^o + ^S^^ (3.25) 



As we noted before, our mechanism for increasing Tc relies on the fact that m^g is determined 
by the time averages {g^^A'^)t and {g'^^Af)t, which are nonzero even if {Af^)t = 0. We also 
find it convenient to set g = 1 by rescaling $ and A^. 



For sufficiently low T//i, we expect that (3.24) admits normalizable solutions. Since we 



are interested in the onset of the instability, we look for a marginally stable mode - namely 
a solution $(r) which is normalizable and everywhere smallj^ This also guarantees that our 
approximation of neglecting backreaction is consistent. The boundary conditions then read 



$(r oo) 



|$(r r_i )| < oo . 



(3.26) 



^In the static case this was studied by [12] . 
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Notice that the regularity condition at the horizon is different from normahzabihty; this 
condition will play an important role shortly. 

The existence of marginally stable solutions can be reduced to a quantum mechanics 
problem. Redefining 



rf{ry/' 



.|)(r) 



sets to zero the friction term, and (3.24) becomes 



Here we have defined 



V{r) = 



eff _|_ 1 _|_ /' 



(3.27) 



(3.28) 



(3.29) 



/ 2 / r/ ^ 

and primes denote derivatives with respect to r. Therefore, we are looking for bound states 
of the Schrodinger equation in the central potential V{r), with energy = 0. Note, however, 
that the normalization condition on $ is distinct from the usual one for the wavefunction 
in quantum mechanics. (It is possible to further redefine variables to obtain the standard 
normalization condition, but we do not do so here for brevity.) 

The strength of the negative term in the potential is controlled by the dimensionless 
parameter 

/i^ + 5/iV2 



A 



(3.30) 



As we increase A we expect to find a critical value Ac at which the first bound state appears 
with E = This configuration corresponds to the onset of the superconducting transition 
and determines T/. 



Tr 



47r 



A. 



1/2 



(3.31) 



where we used (3.3) 



The critical eigenvalue will be a function 



Ac 



/i(A+,(5/iVw^), (3.32) 

which are the remaining parameters appearing in (3.29 )j^ 

We performed a numerical analysis of the first marginally stable solution and the function 
h. The numerical routine imposes the boundary conditions 

$(r+) 



$'(r. 



-A+(A+-3)- 



(3.33) 



Naively one might expect the behavior of (3.28) to depend on four dimensionless parameters: A, 5^ /uj , 



A+, and Tj^/L. However there is an additional scaling symmetry, under which r and r_(_ have charge 1, and 
everything else is neutral, ruling out ?'+ / L. This corresponds to the fact that on the dual field theory side 
Tc should not depend on L. 
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and then scans over increasing values of A until (3.24) admits the first nontrivial solutionj^ 
We find that h is well approximated by a linear function in ^yU^/cu^, so we arrive at the 
expression for the critical temperature 




/^o + ^/^^ 

7o + 7i¥ 



1/2 

(3.34) 



Here 70 and 71 are functions that increase monotonically with A_|_; 70 has been analyzed 
in . An example of this behavior and the fit are shown in Figure |2| for the case l\+ = /. 




Figure 2: Behavior of /i as a function of for = 2. The blue line shows the numerical 

result, and the red line is the best fit, whose equation is also shown in the plot. 



The formula (3.34) for reveals some interesting properties of our system. First, in both 
the static and time-dependent cases Tc decreases as we increase A+. The reason for this is 
that increasing A+ increases the scalar mass, which disfavors an instability. As an example, 
for a fermion bilinear O = iptp, the operator dimension is A+ = 2 (at weak coupling and, 
with suitable supersymmetric protection, also at strong 'tHooft coupling). In this case, the 
critical temperature becomes 

3/ l4-^¥^ Y\ (3.35) 
4vr \^ 16. 75 + 1.48^/ 

We also see that a slowly-varying chemical potential degrades the superconducting prop- 
erties of the system, by lowering the value of T^, however, above a certain critical value 
u) > oOc, the time-dependent source enhances superconductivity. Demanding Tc > Ti°\ 
where Ti°^ = J;:/^o/\/7o is the transition temperature in the static case, we obtain 



cOr. = W — /io = ^V^^i°^ ■ (3-36) 
7o 3 
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This is similar to the static case calculations in [7 . We thank S. HartnoU for discussions on this point. 
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A critical frequency ujc also arises in BCS superconductors, as first recognized by Eliash- 
berg [1] who found 

ujI ~ - (3.37) 

^0 

where tq is the relaxation time of the fermions, and A is the gap. Roughly speaking, one must 
excite the system at a rate fast enough to compete with the dissipation of the quasiparticles. 



4 Analysis of the Condensed Configuration 

Having determined that the instability is enhanced by virtue of time dependence, we now 
consider the solution where the scalar field condenses, for temperatures T < T^. We first 
describe the basic features of the condensed solution and analyze the order parameter {(D) oc 
0+. Next, we check that our classical time- averaged analysis provides a good approximation 
to the behavior of our system. 



4.1 Description of the Condensed Solution 



So far we have determined that T^. has increased according to the classical time-averaged 
solution, but have not exhibited the condensed solution itself for T < T^.. In the full solution, 
the gauge fields are deformed away from (3.22) and (3.19), responding to the nonzero 
in the last term of each of equations (3.9) and (3.10). These effects have to be taken 



into account in order to determine the condensed configuration. Indeed, the backreaction 
5A^ oc generates nonlinear terms in $. Balancing them against the tachyonic mass 
produces a nontrivial profile for $ that we now study. 

In the large frequency regime the deformation of A^. is suppressed by additional powers 
of l/w and may be neglected. The deformation of At is sourced by the term <^^At in the 
Maxwell equation. This is a rapidly oscillating term, controlled by /i(t), but in order to 
capture the backreaction on the time- averaged solution for we must include corrections to 
the rapidly oscillating solution for At. As in our discussion above of the unperturbed solution, 
the interaction of the gauge field with $ is proportional to A\, which does not vanish on 
time-average; therefore, we must solve the Maxwell equation for the back reaction including 
the fast oscillations, and then solve the $ equation in a time-averaged approximation. 

Let us now consider the condensed solution. It is important to analyze the behavior of 
the order parameter near T f» first, as this characterizes the phase transition in the 
dual superconductor. In this limit, the backreaction on the gauge field is small, so we only 
need to keep the leading effect. To work this out, we vary a little bit away from Tc, and the 
resulting scalar field must be small: $ ~ e$0; where $o is the normalized solution to equation 



(3.24) in the non-backreacted case. The leading backreaction to At can be determined from 



its equation of motion as 



A = [/io - 5/isin(u;t)] | - ^) 



h{r-) — h{r^)- 



(4.1) 
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where h{r) is defined as 



hir) 



dr' 

r'2 



oo r, //2 

dr" 



fir") 



ll\2 



%{r") 



(4.2) 



and the h{r^)r^/r term makes sure that At vanishes at the horizon. Here the term of order 
is the deformation of At due to the backreaction. Let us note that this deformation is 
negativep] therefore making a positive contribution to in (3.24) that is of order e^. This 
must be comp ensated by varying T away from Tc by an appropriate (negative) amount, which 
changes (3.24) by terms of order T — Tc. Therefore, e must be proportional to (T — T^)^/^. 



Putting back the dimensional parameters, the value of the order parameter in the large 
frequency limit becomes 



{O) 



T/T, 



,1/2 



(4.3) 



with ca+ a monotonically increasing function of A+. The main features of this result are 
the overall factor r^+, which reflects the dimension of the operator O dual to $, and the 
dependence (1— T/Tc)^/^. This is the expected critical exponent from a mean-field description 
of the superconducting transition. This expression for O is formally the same as in [7J, but 



with Tc now enhanced by time-dependence, Eq. (3.34). 



On the gravity side we can go beyond the LG approximation and study the system at 
T <^ Tc, where it is found that the order parameter has the limiting behavior {O) oc T^+. 
We refer the reader to [7] for more details. 



4.2 Classical and Quantum Corrections to the Time-Averaged Approximation 

We find that the time-averaged approximation is a valid approximation at least for large 
frequency u. To check this we first consider the effects of relaxing the time averaging on the 
classical solution. Then we estimate quantum particle production and decay. 



4-2.1 Classical Solutions 



Let us consider the time-dependent equation of motion (3.11) for the scalar $ in the back- 



ground gauge field (3.21) and (3.22). As argued earlier, this is the sufficient for determining 
Tc, because at T Tc the backreaction of $ on the gauge field is negligible. We may expand 
the full solution in Fourier series (anticipating that the system will only source modes at 
multiples of the frequency co) 



$(r, t) = (f)o{r) + 01 (r) sin(wt) + 



(4.4) 



""^^This can either be seen numerically, or shown by changing coordinates from r to z = i"^/r. Then h is a 
concave function in z vanishing at z — 0. The deformation is h{z) — h{z = \)z and therefore non-positive 
for any z or r. 
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The time-averaged approximation allows us to neglect 0i (and all higher Fourier modes). 
Let us work out the corrections if we had included (pi. At this order of Fourier expansion, 
01 satisfies the following ordinary differential equation 



1 d 

J.2 (J^j. 



dr 



7 



1 - 



0, 



(4.5) 



where m^g is defined by (3.25). This is the same equation as (3.24) save for the addition of 
the last two terms. The last term is a source for 0i; without it we would have the static case 
with vanishing 0i. Therefore, we expect 0i to be linear in this source term (at least when it 
is small). At large frequency u the last two terms compete with each other, giving 



1 



5o < n — 00 



00^ 



00^ \ r 

We can also work out the equation satisfied by the zero mode (po'- 



1 d 

^2 



dr 



2 , , r/^o5/i 



/ 



r+y 

r / 



0, 



(4.6) 



(4.7) 



where the last term is a source. This term is the leading correction to the time-averaged 
approximation. We can estimate its effect by taking its ratio to a typical term in the equation, 
e.g. the /Xq -|- term inside ml^. According to (4.6) this ratio is 



HoSiMpi 



< 



(4i 



So as long as the frequency u is much larger than either g/ig or q6fi, the correction to the 
time-averaged approximation is small. 



4-2.2 Stability and (Un)partide Production 

So far we have considered the behavior of the classical solution in the presence of a time 
dependent chemical potential. Time-dependent perturbations of a field theory generically 
lead to additional quantum effects such as particle production (or in the case of conformal 
field theory, 'unparticle' production) as well. In this section we would like to estimate these 
additional effects of the time-dependent chemical potential Sfi(t) in our system, expanded 
around the classical solution ^^^\r,t), A^^^ for the condensed scalar $ and the gauge fields 
Afj_. To estimate particle production on the gravity side we must study the perturbations 
without time-averaging the chemical potential. 

Let us work in a locally flat frame with coordinates T~tr/L,X~xr/L,i?~L log(r / L) . 
Expanding about the classical solution A^/i'\ we expect the scalar field perturbation 5$ 
to be stable in the time averaged case as in the case of static holographic superconductors. 
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Consider a toy model for our problem encoding this feature, in which the Fourier modes 5$^ 
locally satisfy an equation of motion of the form 



—— + n%{t)5<^ = 0, with nl{t) = —5fi^ir]o/2-cos\uLT/r)) + K''. (4.9) 
al ^ 

Here we take ?7o > 1 to model the fact that the lightest perturbation is stable in the time- 
averaged case. We read off the coefficient L^5/i^/r^ of the perturbation from the action 
expanded in our locally flat frame. 

The first question that this raises is whether a tachyonic mode develops during the times 
when r2^(i(:) < (which occurs for sufficiently small proper momentum K\ Such a tachyon 
grows at most exponentially as ^^i^^'^/^i' = ^Sfit/V2^ rjj^g timescale ~ for this to grow 

appreciably is much slower than the period ~ over which the mode remains tachyonic, 
so the perturbation does not grow exponentially (see e.g. ^'^^ ^ recent discussion of this 
effect). 

The next question is whether there is significant particle production, and if so what its 
consequences are. In particular, the greatest production occurs near the resonance frequency. 
At the Gaussian level for the perturbations of the charged scalar, to calculate particle 



production one would read off the Bogoliubov coefficients from the solutions to (4.9). At 
leading order in the perturbation, the pair production can be computed using ordinary 
quantum mechanics by analyzing the two-state system consisting of the ground state and 
the two-particle state. For energy levels which are separated by an energy difference from 
the vacuum of near the resonance frequency, such a two-state system oscillates between the 
vacuum and the excited state [T3], with a proper frequency of order Usiow = Sfi'^L/ur (see e.g. 
[IB] for a recent discussion of this effect applied to particle production). So for sufficiently 
short timescales this is a small effect and the probability amplitude is of order the nontrivial 
Bogoliubov coefficient whose square gives the number density of produced particles 

The resulting proper energy density in produced 5$ particles near resonance is of order 

ps^ ~ oo^pioOsion^Tf ~ ujI{6^^) (4.10) 

where we defined the proper frequency Up = uL/r. Given this result for (5$^), we can 



plug it into (3.9) to determine on what timescale hackreacUon the produced particles source a 
deformation of At which competes with its original magnitude ~ This timescale depends 
on the radial position in the bulk at which we analyze it. Near ~ T^L^, it turns out to be 

tbackreaction ~ Tc/Sfj'^ ~ l/SfX. (4-11) 

Although parametrically longer than in our setup, this timescale is very short if we 
think of our system as a toy model for a real material. 

However, the quartic coupling leads to decay of 5$ modes into gauge bosons. This 

rate F is of order Up = uL/r, and is much greater than cOgiow above precisely for o;^ ^ Sfi^, 
our regime of interest. So in our regime of parametrically large u'^/Sfi^ the produced 5$ 
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particles decay long before their number builds up to the point where 5$ backreacts on the 
solution. The decay products consist of a superposition of high-frequency electromagnetic 
waves. These themselves can annihilate back into $ particles and also can mix with other 
sectors (for example, neutral moduli fields, other gauge bosons, and other charged fields). 
The net effect of this on w^e// scalar condensate depends on the details of these 

sectors, but it is clear that with sufficient mixing to modes neutral under the gauge symmetry 
corresponding to A^, this backreaction can be small. 

Note that the energy density of produced particles can remain well below the gravity- 
side Planck density, and also well below the energy density M|,;^^^^/L^ which would source 
significant corrections to the background geometry. This means that the system is not 
substantially heated by the process until the energy builds up to an extent that the probe 
approximation breaks down, something we can formally push parametrically as far away 
as we wish. In real systems, it has been argued [1][5] that a steady state can be reached, 
balancing dissipation and time-dependent excitation. It would be interesting to explore this 
upon relaxing the probe approximation in the strongly-coupled systems we are studying. 



5 Discussion 

It is interesting to compare and contrast the stimulated superconductivity in the holographic 
theories of §§3 and 4 to the theory reviewed in §2. In the former, the constituents of the 
scalar condensate, to the extent that that notion applies, are strongly bound. In the latter, 
the kinematics of the Fermi surface plays a key role. In both cases, the effect is sensitive to 
the square of the time-dependent perturbation. In the holographic case, this survives time- 
averaging and contributes to m^g. In the BCS case, the condensate determined by equations 



(2.1 ) and (2.2) is linear in 5n, and 6n is determined by the probability of up-scattering away 
the quasiparticles occupying states near the Fermi surface at finite temperature. This is the 
square of the scattering amplitude induced proportional to the time dependent term in the 
Hamiltonian. Of course, our system is on the face of it very far from BCS theory, manifesting 
only the physics of a bosonic superfluid order parameter. 

In the present analysis, the only limitation on the size of the time-dependent source is the 
probe approximation. This approximation includes arbitrary interactions of the bulk matter 
and bulk gauge fields, which affect each other significantly but do not backreact on the 
background geometry (which is determined by an ambient strongly coupled large-N sector). 
This is a rather different regime of control than that obtained in the previous works such as 



[lj[5j, which were limited to a treatment perturbative in 6n (2.2). One question raised in 
those works was what happens beyond that approximation; in the holographic examples, the 
beneficial effects of time dependence get stronger in this regime, as long as we simultaneously 
increase u, keeping u/Sfj, large. 

Given this, it would also be very interesting to go beyond the probe approximation, 
and study the effects of our time-dependent source including gravitational backreaction. 
The static superconductor beyond the probe limit was analyzed in [20], where it was found 
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that the superconducting instabihty persist also in the backreacted solution. A first step 
in this direction would be to define a consistent time-averaged approximation including 
deformations of the metric. Such an analysis would extend our mechanism to systems that 
do not admit a probe limit, and may reveal new constraints on the increase on by time- 
dependent effects. 

So far we have focused on a sinusoidally oscillating chemical potential. Perhaps it is 
more reahstic to consider irradiating the system with electromagnetic radiation. This can be 
analyzed similarly; it introduces various contributions to m^yj through and (^spatja/)- 
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